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I perform a complete classification of 2d, quasi- Id and Id topological superconductors which orig- 
inate from the suitable combination of inhomogeneous Rashba spin-orbit coupling, magnetism and 
superconductivity. My analysis reveals alternative types of topological superconducting platforms 
for which Majorana fermions are accessible. Specifically, I observe that for quasi-ld systems with 
Rashba spin-orbit coupling and time-reversal violating superconductivity, Majorana fermions can 
in priniciple emerge even in the absence of magnetism. Furthermore, for the classification I also 
consider situations where additional "hidden" symmetries emerge, with a significant impact on the 
topological properties of the system. The latter, generally originate from a combination of space 
group and complex conjugation operations that separately do not leave the Hamiltonian invari- 
ant. Finally, I suggest alternative directions in topological quantum computing for systems with 
additional unitary symmetries. 
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Introduction 

The breakthrough concept of emergent Majorana fermions (MFs) in artificial topological superconducting (TSC) 
devices, pioneered by Fu and Kane motivated a number of recent experiments 0] that have already provided 
the first promising results. The two authors demonstrated that the helical surface states of a three-dimensional 
topological insulator Q, with proximity induced superconducting gap A, behave as a time-reversal (T) invariant TSC. 
When a magnetic field is applied perpendicular to the topological surface, a single MF appears per superconducting 
vortex. Shortly after this proposal, it was recognized that the catalytic presence of spin-momentum locking could be 
alternatively provided by the spin-orbit interaction of a Rashba semiconductor [H, Q . In the case of a semiconducting 
wire [H, fabricated for instance from InSb, a Zeeman energy /i g |B| is sufficient to lead to MFs localized at the edges. 
The existence of confined and protected edge MFs is crucial for applications in topological quantum computing |(| 0] ■ 
A pair of MFs defines a topological qubit, which is in principle free from decoherence and protected against noise, 
in stark contrast to traditional spin Q and superconducting qubits Q . Furthermore, edge MFs can also give rise to 
unique transport signatures [Tol - fl3j . such as the usual [Io|, EH or the magnetically controlled [l2[ 47r-Joscphson effect. 

In the case of a semiconducting wire @ with proximity induced superconductivity, the system transits to the topo- 
logical phase when the criterion fi s \B\ > \J + | A| 2 is satisfied (/i defines the chemical potential). The concomitant 
requirement of a high magnetic field can impede the nanofabrication of the device or restrict the possible manip- 
ulations on the MFs. In fact, several proposals concerning quantum information processes rely on the application 
of strong antiparallcl magnetic fields at the nanoscale level [1J], something not easily realizable in the lab. As an 
answer to these obstacles alternative types of engineered TSCs have been put forward, which support MFs without 
the necessary presence of spin-orbit coupling or the application of a magnetic field |15l — Fl9j | - In most of these proposals, 
some kind of inhomo gen eous magnetic order coexists with intrinsic or proximity induced superconductivity. For some 
of these models [l9[ it has been argued that there is a mapping to the case of the semiconducting wire-based 
TSC mentioned above. Nonetheless, this mapping does not necessarily imply an exact topological equivalence. In 
order to discern if two topological systems are equivalent, the topological class and the role of dimensionality must 
be simultaneously taken into consideration. 

In this manuscript I present a complete topological classification of low-dimensional TSCs that support MFs and 
originate from the combined presence of inhomogeneous Rashba spin-orbit coupling v(r), magnetism M(r) and su- 
perconductivity A(r). My primary goal is to shed light on the topological connection between different existing 
proposals for engineered TSCs and in addition to propose alternative advantageous platforms. For my analysis I will 
consider 2d, quasi-ld and Id systems. The quasi-ld case is obtained from the strict 2d case by the inclusion of a 
confining potential V(r). My study provides new engineered TSC that are experimentally accessible. For example, I 
find that in the case of a quasi-ld system with Rashba spin-orbit coupling and time-reversal symmetry (T) violating 
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superconducting order parameter, MFs can emerge even in the absence of magnetic fields or any type of inhomogc- 
neous magnetism. In addition, I examine the effects of dimensionality on the robustness of MFs through separating 
the systems under investigation into weak and strong engineered TSCs. Furthermore, I illustrate that so far over- 
looked discrete symmetries, that I shall refer to as "hidden" symetries (0), distinguish models previously considered 
as topologically equivalent. Generally, hidden symmetries can be either unitary or anti- unitary and result from a 
combination of space group, time-reversal or other internal symmetry operations that when considered separately do 
not leave the Hamiltonian invariant (e.g. (20l. [2l|). Here I discuss two examples of hidden symmetries: i) a unitary 
hidden symmetry resulting from the combination of a reflection and a translation and ii) an anti-unitary symmetry 
resulting from the combination of time-reversal and translation operations. Finally, I also discuss new topological 
quantum computing (TQC) perspectives that appear when hidden symmetries are present. 

At this point, I give a brief description of how the several sections are organized. In Section I, I provide a short 
introduction to Majorana fermions and introduce the general Hamiltonian that describes the systems of interest. In 
Section II, I shortly review the topological classification methods with special focus on the situations where additional 
unitary and anti-unitary symmetries are present. In Section III, I present an overview of my main results (Table [TTJ) 
concerning the classification of TSCs when all possible spatial symmetries are broken. I further discuss how the 
emergence of hidden symmetries can modify Table [TTJ In Section IV, I provide a detailed analysis and justification 
of the results presented in Section III. In Section V, I present two specific examples of systems characterized by a 
hidden symmetry and study the impact of the latter on the topological properties. In Section VI, I discuss how the 
presence of hidden symmetries can be useful for developing topological quantum computing protocols and suggest 
possible candidate systems that could be used for this purpose. Finally, Section VII summarizes my main results and 
related conclusions. 



I. MAJORANA FERMIONS AND MODEL HAMILTONIAN 



In condensed matter physics MFs are not fundamental particles [22J but excitations of a many-body system [23, [24| . 
Essentially, what we define as MFs are the operators 7 Q (a is just a label) which satisfy {-f a , -fp} = S a ^l (I the identity 
operator) and constitute zero energy eigen-operators of the Bogoliubov - de Gennes (BdG) Hamiltonian. Since MFs 
are hermitian they can be described by the following general expression 



dr 



^, a (r)ih(r) + ul ta (r)ipi(r) + u t , Q (r)^|(r) + u ita {r)tj}[{r) 



(1) 



where ipt ( r ) /Vv ( r ) correspond to the creation/annihilation operators of an electron in 2d coordinate space r = (x, y) 
with spin projection a =j" , J,. Notice that MFs require linear combinations of electronic operators and their hermitian 
conjugates. This limits us to systems in which the spin-quantization axis is fixed. Notice that for a system with 
spin-rotational symmetry, the application of a homogeneous magnetic field breaks the latter symmetry but the spin- 
quantization axis can be always redefined. In this case, MFs are not accessible. As a matter of fact, MFs can appear 
only in systems with spatially inhomogeneous magnetism, spin-orbit coupling or spin-triplet superconductivity. 

In this work I focus on systems that satisfy the above requirements and arc cither microscopically or phenomeno- 
logically described by the following bulk Hamiltonian 

— — fj, + V(r) M(r) ■ a + M^M" 
2m 2 



H 



dr ijj\r) 



dr 



(2) 



where ffi(r) = (V'j( T *) V'H'*))) ar e the spin Pauli matrices, v(r) is the spatially dependent strength of the Rashba 
spin-orbit coupling, M(r) corresponds to a magnetic-field or a magnetization profile and A(r) defines a spatially 
varying superconducting order parameter. Notice that in some sense the above Hamiltonian is overcomplete, since it 
covers all the cases that we will consider, without implying that all the terms are simultaneously required for obtaining 
a TSC. Furthermore, at the level of my topological classification, the origin of the involved terms is unimportant. 
However, I have to remark that when I will discuss specific cases I will concentrate on engineered TSC, which for 
instance involve conventional types of magnetism and mainly proximity induced superconductivity (25j . This means 
that I will not consider here the cases of unconventional (2(| density waves [27l - [3l| or superconductors [Hj], although 
some of our conclusions could also apply to these systems. 

Since for the situations considered in the present study the spin-quantization is always fixed, I will employ the 
following spinor 



$ f (r) = U\{r) ,^l(r) ,^ t {r) ,-ifj^r) 



(3) 
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and use the r Pauli matrices in order to represent matrices in the Nambu particle-hole space. With the introduction 
of the above enlarged spinor the Hamiltonian can be rewritten in the following compact way 

^i = \ J dr&(r)H(p,r)V(r), (4) 

where H(p,r) corresponds to the BdG Hamiltonian. Notice that the factor of 1/2 is crucial for avoiding double 
counting of the degrees of freedom, since the above spinor does not obey to the usual fermionic commutation relations. 

II. TOPOLOGICAL CLASSIFICATION PRINCIPLES 

Before discussing the possible topological phases arising from our model Hamiltonian, I will briefly review the basics of 
how to classify topological systems. My goal is to first highlight a key point which is crucial for classifying TSCs and 
then demonstrate how this can provide further topological insight concerning previously studied s yste ms fig . Il9j . This 
key point is that topological classification of systems following the recently developed methods [33ll34^ , is conducted for 
irreducible Hamiltonians, for which one cannot find any unitary operator O u satisfying [H(p, r), O u ] = 0. If there is a 
number of these type of operators, we can block diagonalize the Hamiltonian and topologically classify each sub-block. 
Of course, this is not the only route to study topological properties, since one can also directly construct topological 
invariants for reducible Hamiltonians [35j | . Nevertheless, studying irreducible Hamiltonians provides a transparent 
analysis of the topological classes. Here I am interested in "hidden" unitary discrete symmetry operators satisfying 
the cyclic property C™ = I, with n 6 Z. In the simplest case n = 2, we can block diagonalize the Hamiltonian into 
two sub-blocks labelled by the eigenvalues ±1 of O u , leading to a direct sum of the form 7i + (p, r) © (p, r). Notice 
that because of the discrete symmetry O u , both sub-blocks are constrained to belong to the same symmetry class. 
However, the two sub-systems do not necessarily reside in the same topological class. 

The symmetry class and the related accessible topological phases of an irreducible Hamiltonian are defined by the 
possible presence of three specific types of discrete symmetries. The first two correspond to a generalized time-reversal 
symmetry effected by the anti-unitary operator and a charge conjugation symmetry effected by an anti-unitary 
operator 3. If is a symmetry of the Hamiltonian, it satisfies [H(p, r), 0] = Q~ lr H(p, r)0 = +H(p, r) while in 
the case of charge-conjugation we instead have {H{p, r),E} = => 5~ 1- H(p, r)3 = —W(p,r). If and 3 constitute 
symmetries of the Hamiltonian at the same time, then the Hamiltonian additionaly satisfies {H(p, r), 03} = where 
the combined 03 operator is unitary and is termed chiral symmetry operator n. The inclusion of n completes the set 
of symmetries that are required for determining the symmetry class of an irreducible Hamiltonian. In fact, in order 
to cover all possible symmetry classes, we have to take into account the case in which a unitary chiral symmetry may 
exist without the necessary presence of and 5 symmetries. 

In this work, I also provide an example of an anti-unitary hidden symmetry O a . In this case O a constitutes an 
additional generalized time-reversal symmetry which modifies the initial symmetry class of the system, instead of 
splitting the latter in a direct sum of identical symmetry classes as for the unitary analog O u . For instance, if a 
system is initially in class D, then the emergence of an anti-unitary hidden symmetry O a with 0\ = +1 will change 
its symmetry class to BDI. 

For the cases under consideration, the BdG Hamiltonian enjoys a charge-conjugation symmetry 3 = t x IC, where K. 
defines complex conjugation. Since S 2 = +1, we obtain only the following three allowed symmetry classes presented 
in TableQ] BDI, D, Dili or their direct sums BDI©BDI, D©D, DIIIffiDIII in the presence of a hidden symmetry O u . 
Notice that the classes BDI and Dili are characterized by a time-reversal symmetry with 2 = +1 and 2 = —I, 
respectively. In the first case, symmetry implies that the Hamiltonian is real while in the second that there exist 
a Kramers-type degeneracy leading to doublets of solutions. Below I examine the minimal cases that can lead to a 
symmetry class supporting MFs. For completeness I will also shortly discuss previously studied models. 



TABLE I: Symmetry classes of topological superconductors supporting Majorana fermions. For O and H, ±1 corresponds to 
the result of O 2 and E 2 . For II, I denotes that the symmetry is present. Conversely, implies that the corresponding symmetry 
is broken. 
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FIG. 1: Point group symmetries for quasi- Id and strictly Id geometries of a topological superconductor (depicted with blue). 
In the quasi-ld case, inversion symmetry along the z— axis, ah, is broken. 

III. RESULTS: POSSIBLE PHASES OF ENGINEERED TOPOLOGICAL SUPERCONDUCTORS 

In the present section I carry out a thorough analysis of the accessible TSC phases that follow from the Hamiltonian 
of Eq. ([2]). For the strict 2d and Id cases I will consider that V(r) = 0. To analyze the quasi-ld case, I will always 
assume the presence of a confining potential V(y). For topological computation applications based on edge MFs the 
quasi-ld and pure Id setups are mostly relevant. The possible unitary symmetries that can appear for these systems 
originate from the point group Q and translation operations t a with t a r = r + a. Let me now focus on the point group 
symmetries for the quasi-ld and pure Id geometries, which I also graphically demonstrate in Fig. [T] The point group 
for a quasi-ld system confined in the xy-plane is C^. This symmetry group includes a C2 7r-rotation about the z-axis 
(r — > — r, z — > z) and two o~ v reflection operations <j xz (y — > — y) and a yz (x — > —x), where the indices correspond to 
the mirroring plane. Notice that the reflection symmetry operation er^ = <j xy (z — > —z) is broken in C2,,. In the strict 
Id case we are left only with a yz . For random v(r), M{r) and A(r) all the aforementioned symmetries are broken. 
Nevertheless, for special spatial profiles of the latter functions, a hidden symmetry can emerge, which consists of these 
basic symmetry operations or other already broken symmetries such as ov 

In Table [HI I present the topological classification for the Hamiltonian of Eq. ([2]) where all possible unitary symme- 
tries are broken due to the spatial dependence of v(r), M(r) and A(r). I demand that |Af(r)| ^ and |A(r)| ^ so 
as to avoid any gap closings that could lead to a macroscopic coexistence of different topological phases throughout 
the volume of the material. I also have to remark that in the case of a translationally invariant system, we can 
transfer to fc— space in order to calculate topological invariants. If transnational symmetry is broken, then analysis of 
the topological properties in coordinate or momentum space exhibits the same complexity. Of course, there can be 
also cases where topological properties in combined (r, fc)-space can be relevant [35l. |36|. 

One of the most important results is Case II in Table [Til where the simultaneous presence of Rashba spin-orbit 



TABLE II: Accessible topological superconducting phases due to the combined presence of inhomogeneous Rashba spin-orbit 
coupling v(r), magnetization M(r) and superconducting order parameter A(r). The presence (absence) of the aforementioned 
terms is indicated with + (-). The resulting symmetry class depends on the behaviour of the magnetic and superconducting 
Hamiltonian terms under complex conjugation /C, since the Rashba spin-orbit coupling term always preserves T. For a term 
that is already present (+), we denote the case of preserved (broken) complex conjugation as K, — I (0). Notice that symmetry 
classes that lack a strong topological invariant for the corresponding dimensionality are shown with italics and correspond 
to weak topological superconductors. In the presence of a unitary hidden symmetry O u with the cyclic property 0™ = I, 
a symmetry class TC splits into n identical sub-classes ©„TC. However, the sub-systems do not have to reside in the same 
topological class. 
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coupling and inhomogcneous superconductivity can lead to MFs in a quasi-ld system, without the requirement of a 
magnetic field. In fact, an r-dependent superconducting phase originating from a supercurrent falls into this case, 
constituting an experimentally prominent route towards MFs. As far as the table is concerned, the possible phases 
are essentially classified by the behaviour of the magnetic and superconducting Hamiltonian terms under K,. 



IV. ANALYSIS OF THE POSSIBLE TOPOLOGICAL PHASES IN THE ABSENCE OF UNITARY 

SYMMETRIES 



In this section I provide the detailed topological classification for the cases presented in the Table [Til Notice that for 
the present discussion the spatial dependence of the terms involved is considered random, unless explicitely stated. 



• Cases I and II 

In the following paragraph I will focus on the Cases I and II that are characterized by the presence of inhomo- 
geneous Rashba spin-orbit coupling v(r) and superconducting order parameter A(r). The TSCs belonging to 
these cases are described by the following Hamiltonian 



2m 



M + V(r) 



{v{r),P x T z (Jy ~Py<T x } . , s A / \ 

^ y - A R (r)TyCTy- A^(r)T x <j y . 



(5) 



The Rashba spin-orbit coupling term is odd under inversion symmetry along the z-axis ah-, while it is even under 
the usual time-reversal symmetry T. If the superconducting term is also invariant under T or equivalently )C, 
since we are dealing with a scalar superconducting order parameter, then A(r) = Asr(j") and the full Hamiltonian 
is characterized by the generalized time-reversal symmetry 6 = ia y JC that coincides with T. 

2d system: In the 2d case, the particular system belongs to the symmetry class DHL Since 9 satisfies 2 = —I, 
with I the identity operator, we expect boundary MF Kramers doublets. Class Dili possesses a strong Z2 
topological invariant in 2d. The presence of also leads to a chiral symmetry with n = r x a y . In the case where 
the superconducting order parameter has an additional imaginary component, T is broken and the system 
transits to class D. Class D has a strong Z invariant in 2d and consequently this system constitutes a strong 
TSC in both cases. 

In order to analyze the symmetry properties in a more transparent manner, I will consider without any loss of 
generality, the following form for the superconducting order parameter A(r) = Ae lJ r . The particular profile, 
constitutes the simplest representative of T violating superconductivity and can be viewed either as the result 
of the spontaneous formation of a Fulde-Ferrell [13] phase with modulation wave- vector J or the consequence 
of the application of a supercurrent J. The Fulde-Ferrell phase is a special case of pair density waves [38[ that 
have been also recently considered as potential TSCs leading to MFs. On the other hand, the application of 
supercurrents was previously disced in Refs. @. In the latter implementations a supercurrent was viewed 
as an extra knob for tuning the topological phase diagram. 

At this point we proceed with gauging away the superconducting phase ip(r) = J ■ r via the minimal coupling 
p — > p + hVip{r)r z /2 = p + hJr z /2, leading to 

#(P) = ^J-Pl + 
2m 

- Ar a cr y . (6) 



P 

2m 



(hJ/2)' 

2m 



ft + V{r) 



, {v(r),p x T z a y -p y a x } v(r)H 
t z H 1 — [J x a y 



JyT z <Jx) 



It is straightforward to confirm that for J = the system belongs to class Dili because T is preserved while 
for finite J the system lies in class D. 

quasi-ld system: In order to investigate the quasi-ld and Id cases I set v(r) = v(x). Furthermore for the 
quasi-ld case I additionaly switch on a confining potential V(r) = V(y). The presence of the confining potential 
lowers the symmetry of the system, permitting anisotropic coefficients for the Rashba terms p x r z a y and p y u x , 
instead of a common v{x). For my analysis I will keep the coefficients equal since the only crucial requirement 
for my study is that they are both non-zero. To achieve confinement, I consider the case of a harmonic potential 
V(y) — muj 2 y 2 /2. This term is translationally invariant along the x-direction and even under C2, o xz and a yz . 
Another option for the confining potential is the infinite wall potential V(\y\ > L y ) = +00. For the choice of 
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the harmonic confining potential, the Hamiltonian reads 

5t\ 



H(p x ,a, at) = —[JxVx + Jy 
2m \ 



mum a — a' 



{v{x),p x } 



T z (7 y - V(x) 



Pi (HJ/2) 2 fc / , 1 



muih a — at 



v(x)h 



(7) 



where I introduced the quantum harmonic oscillator's bosonic creation (annihilation) operator at (a). By 
introducing the eigenfunctions \n) of the number operator N = at a, I obtain the matrix Hamiltonian 



n(p x ) 



2m 



(J X P X I + Jyfr^y) 



2m 



{v(x),p x \ v(x)h . 
[i\ z ) r 2 H T z a y - v(x)h\ y cr x H — {J x a y - J y T z a x ) 



ATyCTy , 



(8) 



that is defined in spin, Nambu and N spaces with 
(n\X J z \s) = 6. 



W2) 2 h ( 1 
/! fa | n 



2m 



,mcj y/n + 16 n>a -i - y/nS n s+1 
jp, and (n\ X y \s) = x j — . (9) 



Since the form of the Hamiltonian is identical to the 2d case and K. X y IC = —X y (similarly to p y ), the quasi-ld 
model also belongs to the Dili class for J = and to class D for J ^ 0. 

Id system: For studying the strictly Id system, I apply the dimensional reduction method to the 2d model of 
Eq. ((6|) and set p y = J y = 0, that yields 



Wipx) = 7— JxPxI ■ 
Am 



Pi , (hJ*/2Y 



2m 



2m 



{v(x),p x } v(x) 
T z a y H — n.J x a„ — H\t„o. 



>x<Jy 



yUy 



(10) 



We observe that for this model we retain our freedom to redefine the spin-quantization axis and as a result the 
above Hamiltonian does not support MFs. 

• Cases III and IV 

In this section I consider TSC phases that do not involve spin-orbit coupling. This implies that at least two 
components of an inhomogeneous magnetization field must be present in order to lock the spin-quantization 
axis, since the latter constitutes a prerequisite for obtaining MFs. For this kind of systems, the Hamiltonian 
reads 



n(jp,r) 



2m 



V(r) 



t z - M(r) ■ (t z ct x , a y , r z a z ) - Asft(r)T y a y - A a (r)r x cr y . 



(11) 



For the specific type of TSCs, the magnetization field M(r) is odd under the usual time-reversal symmetry T. 
However, its behavior under complex conjugation JC is not fixed. If M y (r) = then M(r) preserves /C. This 
leads to the following two possiblities depending also on the behaviour of the superconducting order parameter 
under 1C. In the first possibility the magnetic and superconducting terms are simultaneously invariant under JC 
and a generalized time- reversal symmetry appears with = K accompanied by a chiral symmetry H = t x . 

2d system: In 2d, the system belongs to the BDI class that however is not characterized by a strong topological 
invariant for this dimensionality. Consequently, the specific system corresponds to a weak TSC, since under 
special circumstances one can generally define weak invariants. The second possibility involves the breaking 
of JC by cither one of the terms. In the latter case, the Hamiltonian belongs to class D which has a strong Z 
topological invariant in 2d. 

quasi-ld system: For the particular study I will consider for convenience that A(r) = Ae lJ r . As previously, 
I gauge away the superconducting phase and obtain the equivalent model 

t z - M(r) ■ (t z <t x ,<t v ,t z <7 z ) - ATyOy . (12) 



n'(p) 



2m 



J pi 



P 2 , (W/2) 2 
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For effecting confinement I will employ once again a harmonic oscillator's potential V(y) — muj 2 y 2 /2 and we 
also have M(r) = M(x, a + a^). Following the same steps as in Cases I and II, I obtain the Hamiltonian 



U'{f>x) = ( JxPxI + JyhXy) + - n\ J z ) t z - M(x) ■ {r z a x ,a y ,T z a z ) - A Ty a y , (13) 



where M(x) is a symmetric matrix defined in \n) space. If J = and M y (x) = 0, K. is a symmetry of the 
Hamiltonian and the system belongs to class BDI. Instead, if M y (x) ^ 0, the system belongs to class D. For the 
special case where M(r) does not depend on the y-coordinatc, i.e. M(r) — M(x), M{r) becomes diagonal 
and can be divided into an infinite number of sub-spaces labelled by n yielding 



Pi fc ( 1 
u + nuj \ n H — 

2m p V 2 



t z - M(x) ■ {r z a x ,a y ,T z cr z ) - Ar y a y , 



(14) 



which leads to the total symmetry class ©„BDI. If M(x) violates JC we obtain a direct sum ©„D. By allowing 
a finite J we also violate JC. Specifically, if J = ( J x , 0), that corresponds to the case A(r) = A(x), the system 
resides in the class © ra D. However, if J = (0, J y ) the system belongs to class D, due to the simultaneous presence 
of A y and X z in the Hamiltonian of Eq. (fT3]) . that does not allow the decomposition in n-sectors. In Table [TT1 
the general case where M and A depend on both coordinates is depicted. 

Id system: By dimensional reduction on the Hamiltonian of Eq. (TT21) we obtain the following pure Id model 

t z - M(x) ■ (T z a x ,a y ,T z a z ) - AT y a y . (15) 



2m 



pi , (hJ x /2y 



2m 



2m 



If M y (x) = and J x — 0, JC is conserved and the system belongs to class BDI. Instead, if one of the previous 
terms is non-zero, the Hamiltonian is not real any more and it falls into symmetry class D [I?}. 

• Cases V and VI 

In the last part of this section I complete the possible cases by considering the situation where all the terms of 
Eq. ([2]) are present. The latter equation in combined Nambu and spin spaces reads 



H(p,r) 



Pi 

2m 



- A* + V(r) 



{v(r),p, 



jCT x } 



- M(r) ■ (r z a z ,a y ,T z a z ) - A^(r)T y (jy 



A^(r)T x a y . 



(16) 



When magnetism and Rashba spin-orbit coupling coexist, the accessible topological phases constitute an overlap 
of the previously examined separate cases. Therefore here we will investigate what are the consequences of 
the addition of magnetism in Cases I and II for different dimensionalities. Earlier, we observed that when 
magnetism is not present, there are two possible scenarios depending on the behaviour of the superconducting 
order parameter under JC. 

2d and quasi-ld systems: For the specific cases, if A(r) = Asr(t*) the system resides in the symmetry class 
Dili being invariant under T . If M{r) is introduced, T will be broken and the system will transit to class D. 
If A(r) is complex, the system is already in class D, and consequently the inclusion of magnetism leads to no 
additional effects. 

Id system: For pure Id systems the presence of a magnetic order is crucial and leads to new TSC phases. The 
Id descendant of the above Hamiltonian reads 

^ / p2 \ -fv(x) Px\ 

H(p x ,x) = ( _ P- J T * H 2 Tz<Jy ~ M ( x ^ ' ( r « £r «i°i/i 7 '* *) ~ Au(x)T y (T y - A^(x)T x a y . (17) 

From Table 1111 wc immediately observe that there are no TSC phases in the absence of magnetism. As mentioned 
earlier, the reason is that the presence of the spin-orbit coupling term T z a y alone, cannot fix the spin-quantization 
axis. Nevertheless, the addition of a perpendicular magnetization field remedies this problem and can lead to 
TSC phases. If M(x) is invariant under /C, i.e. M y (x) = 0, the Hamiltonian is characterized by a generalized 
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time-reversal symmetry = K, and a chiral symmetry IT — t x which permits an integer number of MFs per 
edge [4lJ . The translationally invariant version of this model 

U{p x ) = - [J^J T z + Vp x T z (Ty - M ■ {T z <7 z ,<Ty,T z CT z ) - kTyOy , (18) 

corresponds to the celebrated MF-wire proposal [j| which currently under intense experimental investigation Q 
and concerns a Rashba semiconducting wire in the presence of a Zeeman field and proximity induced supercon- 
ductivity The systems transits to the topologically non-trivial phase when the criterion 



|M| > vV + A 2 , (19) 

is satisfied. Finally, if M y (x) ^ then JC is broken and the system belongs to class D with a I2 invariant 
allowing for a single MF per edge. 



V. EXAMPLES OF TOPOLOGICAL PHASES WITH HIDDEN SYMMETRIES 



In this section I will demonstrate two examples where unitary or anti-unitary hidden symmetries occur for some of the 
TSC phases presented in Table [TT] and demonstrate what are the concomitant modifications of the initial symmetry 
class. 

• Cases I and II in the presence of a single unitary hidden symmetry O u 

Let us now investigate the consequences of the emergence of a "hidden" symmetry due to the special form of 
the Rashba spin-orbit coupling term. As a case study I will focus on the topological properties of the following 
quasi-ld Hamiltonian (in Id we do not obtain a TSC), introduced in Eq. © of Case I 



h 1 T - t , 7 t\ \ , ( Pi ,j\ , {v(x),p x } . ... , v(x)H 



H(px) = 7T~ {JxPxI + Jyfr^y) + I lf~ ~ (J>K T z H T z (J y - v(x)h\ y (J x H ( J x <7 y - J, 



2m v y y ' \2m 

Ar y a y . (20) 



Here we will restrict to the special situation where t^/Qv{x) = v(x + tt/Q) = —v(x). We may readily observe 
in which manner this property leads to an emergent unitary symmetry. The terms of the Hamiltonian that do 
not contain v(x) are invariant under arbitrary translations and under the action of er/j, which in our formalism 
is represented as = ir z a z in spin-space. Since all the terms with coefficient v(x) are odd under ah, the full 
Hamiltonian is invariant under the action of O u = a^t^/Q. The appearance of a hidden symmetry O u leads to 

an additional generalized time-reversal symmetry = O u T and a concomitant chiral symmetry n = O u Tl. The 
emergence of O u modifies the symmetry class of the system by splitting the symmetry classes Dili and D found 
earlier, into DIII0DIII and D0D, respectively. Notice that the possibility of point group symmetry protected 
phases has been recently discussed J42j. A simple example for v(x+n/Q) = —v(x) is v(x) = 2vq cos(Qx + 0) with 
Q = 2q. Here, 9 is considered pinned to a constant value. The modulated spin-orbit coupling term can be viewed 
as an unconventional spin triplet density wave (27j . similar to the Rashba spin-orbit density wave proposed in 
[43j as a potential candidate for the so called "hidden order" , which appears in the non-superconducting regime 
of the phase diagram of the heavy fermion compound URii2Si2. 

The simultaneous presence of the momentum operator p x and coordinate x does not allow for a direct and 
transparent inspection of further topological properties of the system. Nonetheless, it is also possible in principle 
to obtain through some kind of "deformation" procedure (in the topological sense) a model defined solely in 
momentum space that shares the same symmetries and topological properties with the original model. In 
order for this mapping to be meaningful and offer a direct computation of topological invariants, translational 
symmetry must be somehow restored. The presence of a periodic v(x) = 2vq cos(Qx + 9) term, leads to 
the formation of a band structure with a Brillouin zone of length Q since t2^/Qv(x) = v(x). The property 
t n /Qv(x) = — v(x) gives rise to a sublattice structure that will eventually lead to the two sub-blocks of the 
Hamiltonian that become relevant in the presence of O u . Since we are not interested in the full band structure, 
but mainly in restoring translational symmetry, we may expand the field operator in the following fashion 



(21) 
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where ip± q (x) arc slowly varying fields leading to the Hamiltonian 



pi (MY xJ 

2m Im 



VPxT z p z 



+VqPx (cosQT z pxOy - s'm6T z py(Ty) - VqH (cos 9p x X y a x - sinOpyXycr x ) 



-J x (cos Qp x Oy — smOpyUy 



-J y (cos9r z p x a x - sm9T z PyCr x ) - Ar y a y . 



(22) 



with v = hq/m. The above Hamiltonian acts on the enlarged spinor 

%\(x) = (i'l qt (x) , ip\. qi {x) , ^_ qt (x) , ijjt_ qi (x) , ip_ qt (x) , i>_ ql (x) , ip +gt (x) , il) +qi (x)J , (23) 

with the p Pauli matrices acting on ±q space. Notice that terms with p x or p y carry momentum Q. By expanding 
the field operator in this manner, we managed to end up with a translationally invariant Hamiltonian. This 
approximation allows us to readily study topological aspects in momentum space which for the specific case is 
an easier task compared to the required analysis in coordinate space. Nonetheless, it is not a priori ensured that 
the coordinate and momentum pictures are topologically equivalent. If they do, this approximation constitutes 
a suitable deformation procedure for mapping x to k x space topology. 

In order to confirm if these systems belong to the same symmmetry class, we have to study the emerging 
symmetries for the latter model. In this basis the expression for the generalized time-reversal symmetry operator 
simplifies to Q = ia y JC = ip x a y lC', where YJ is a complex conjugation operator not acting on Q or q. The 
presence of p x effects complex conjugation operation K, in g-space, since q — > —q is equivalent to K, = p x lC' . 
Furthermore, within the specific framework t^/q = —ip z , up to an unimportant fc-dependent phase. Notice that 
(—ipz) 2 = —I, i.e. similar to the behaviour of rotation operators for a spin-1/2. This is a direct consequence of 
the fact that the spinor contains the wave-vector q which is half of the wave- vector Q = 2q. We may directly 
confirm that \H q (j> x ), <Jht-n/o\ = [HqiPx), T z pz&z] = 0. When J = 0, the Hamiltonian is invariant under O 
and the presence of O u = Uht^/Q leads to the additional time-reversal symmetry = O u Q = iT z p y a x KJ . The 

emerging chiral symmetries for this model read n = r x a y and n = r y p z a x . As expected, in this case the system 
belongs to the symmetry class DIII©DIII. Furthermore, if J ^ then O is broken and the system transits to 
the symmetry class D©D. To explicitely demonstrate the sub-block structure of the Hamiltonian and the direct 
sum of symmetry classes, I effect the unitary transformation 



A = 



I + Pz^z 



I - PzOz 



(24) 



which transforms the Hamiltonian as follows 'H q ot (p x ) = A^T-L q (p x )A. This particular unitary transformation 
block diagonalizes the matrix r z p z a z , representing the hidden symmetry operation O u , into A^T z p z a z A = t z . 
The transformed Hamiltonian is block diagonal and can be labelled by the eigenvalues of t z , t = ±1, yielding 



U r q 0t T {p x ) = —(JxPx + JyhXy)I+-vJ xPz +T 



f + ^-f^z 

2m 2m 



p z (Tz + rvp x a z 



2m x " 2 

+Vq [jp x cos 6 + hX y sin 9) p y a x + Vq (rp x sin 9 — HX y cos 9) p x a x 

VqH 

H — — [(Jx cos 9 + rJ y sin 9) p x a y - (J x sin 9 - tJ v cos 9) p y a y ] + rAp z a x 



(25) 



I have to remark that the above topological classification conclusions hold for a bulk system. In order to observe 
the two sets of edge MFs one must introduce boundaries that preserve the O u symmetry. Notice however that 
in the present example the preservation of O u requires a specific behaviour under the translation operation 
t-x/Q, which will be unavoidably broken for a finite system. Nonetheless, the case discussed here does not only 
constitute an example of mere academic interest. Although the presence of protected boundary modes [3, H] 
is considered to be the hallmark of topologically non-trivial phases, it does not constitute the unique route 
for diagnosing topological order. In fact, fingerprints of topological non-trivial phases can be also found in 
manifestations of the bulk system. One example is the polar Kerr effect [1^] that characterizes class D chiral 
p-wave superconductors. This experiment can provide a direct evidence of topological order by solely probing 




FIG. 2: Engineered topological superconductor consisting of a metal with spin-spiral magnetic order M(x) = 
(2A'I® cos Qx, 0, 2M® sin Qx) placed on top of a bulk conventional superconductor. The magnetic order is invariant under 
the combined action of time-reversal symmetry T followed by a n/Q translation operation, t n /Q. The presence of the anti- 
unitary hidden symmetry O a = t-w/QT leads to the class BDI©BDI, compared to class BDI when it is absent. 



the bulk response. Similar chiral phenomena emerge in non-superconducting systems. In the latter, apart from 
a similar polar Kerr effect [29|, Etjj . an anomalous thermoelectric Nernst response [3(1 E3| and a topological 
Mcissner effect [3l| constitute additional smoking gun signatures of quantum anomalous Hall phases (class A). 
In fact, topological response survives also in finite temperatures, though exhibiting no quantization phenomena. 
Evenmore, the bulk magnetic response |48| of a quantum spin Hall insulator (class All) can provide alternative 
routes for confirming the transition to the topologically non-trivial phase. 

• Cases III and IV with a single anti-unitary hidden symmetry O a 

Here I will investigate the consequences of the emergence of an anti-unitary hidden symmetry O a on the sym- 
metry class of the Id model of Eq. (fT5|) that was obtained for Cases III and IV. In this way I will be in a position 
to make a connection to previous studies [lH, [l9[. The Hamiltonian of Eq. (|T5j) for J x = 0, has the following 
form 

H(p x ,x) = (J^- - fij t z - M(x) ■ {T z cr x ,<jy, r z a z ) - ATydy . (26) 

First I will focus on the Case III, which belongs to class BDI with O = K., if M(x) is random and M y {x) = 0. 
I demonstrate that the topological properties of the system change with the emergence of a hidden symmetry. 
I now assume that t 7! /QM(x) = M(x + ir/Q) = ~M(x). For this special case, the Hamiltonian is invariant 
under the anti-unitary hidden symmetry O a — t^/gT. Since t^/g is unitary and T anti-unitary, the particular 

symmetry constitutes an additional generalized time- reversal symmetry = O a - In order to gain more insight, 
I will consider the simple spin-spiral ma gne tization profile M x {x) = cosQx and M z {x) = 2M® svaQx 

which is depicted in Fig. [21 Prior studies [l8l [l9j have focused on the case = M®. By expanding the field 
operator as in Eq. (|21[) we obtain 

T z + Vp x T z p z - M®T z p x (T x + M®T z p y a z - ATyCTy . (27) 

Within this framework we have = K = p x K! and = O a = t^/qT = —ip z ia y IC = —ip z ia y p x K! = ip y a y JC' . 
Remember that K! does not act on q and Q. We readily observe that for both generalized time- reversal 
symmetries we have 2 = 2 = +1 leading to the symmetry class BDI®BDI. Note that in Refs. [3 [l9[ it 
was shown that, up to a spatially dependent unitary transformation, the model of Eq. (|2"fi")) is equivalent to the 
MF-wire model of Eq. (fl8]) with the latter belonging to the symmetry class BDI. However, as we showed here the 
particular system belongs to class BDI®BDI. Generally a unitary equivalence does not necessarily constitute a 
topological equivalence. Nonetheless, it is straighforward to demonstrate that the two pictures agree with each 
other. By performing the transformation T-L rot (p x ) = K' { T-L rot {p x )K with A = (p z <J z + <J y ){p y + Pz)/2 we obtain 

t z + vp x T z p y + M®T z p z - M®T z p z a z - ATyp y a z . (28) 



n(p x 



2m 



(hqy 

2m 



rl ro \p x ) 



pj { {hqy 

2m 2m 
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As anticipated, the above Hamiltonian is block diagonal and if we introduce the eigenstates a = ±1 of a z in the 
rotated frame, we obtain 

t z + vp x T zPy + (A/f - <rMf) t zPz - aATyPy . (29) 

Remarkably, each of the above block Hamiltonians is identical to the MF-wire Hamiltonian of Eq. (|18|) with 
effective chemical potential \i — (hq) 2 /2m, spin-orbit coupling strength v, Zeeman field M® — crM® and su- 
perconducting order parameter a A. Each of the blocks will be in the topologically non-trivial phase when the 
criterion 

| MQ - aMf | > y / [ M -(%) 2 /2m] 2 + A 2 , (30) 

is satisfied in complete analogy to Eq. (Q~9j). We observe that the two sub-systems are not necessarily in the 
topologically non-trivial phase, at the same time. Evenmore, if we consider = as in prior studies 
fla . [l9|, only one of the sub-systems can be in the topologically non-trivial phase. In this case, the system 
effectively behaves as a class BDI TSC and this is in accordance with the previous analytical findings. Notice 
that for a bounded system the hidden symmetry will be broken and the multiple edge MFs are not expected to 
be observed even when M® ^ . 



K ot (p*) 



Pi , (HY 



2m 2m 



VI. NEW TQC PERSPECTIVES IN A TSC WITH UNITARY HIDDEN DISCRETE SYMMETRIES 



For the cases that we considered in this work, hidden symmetry involved a specific behaviour under translations. 
However, this type of symmetry will be explicitely broken when boundaries arc introduced. Nevertheless, one could 
look for alternative, robust and tunable, hidden symmetries that arc related to some internal degree of freedom such 
as a valley, orbital or band index. 

Let us now discuss new routes that open up for topological quantum computing when we consider the additional 
presence of a hidden unitary discrete symmetry O u . Generally, two edge MFs 7 a and % combine into a zero-energy 
fermionic mode d = + i7b)/v2 that leads to a doubly degenerate ground state |1) and |0). These two states 
correspond to many-body ground states where the zero-mode is occupied or not, res pec tively. The Hilbcrt space 
spanned by these two degenerate states defines a topological qubit which is in principle [49| protected by decoherence 
due to the non-local binding of the MFs. The fundamental and in fact the only allowed topologically protected single 
qubit operation which we may perform is called braiding 0, 0] and corresponds to exchanging the two MFs in real 
space Fig. [3] After braiding is effected the states |1) and |0) become e 47r / 4 |1) and e~ 47r / 4 |0), picking up a relative 7r/2 
phase. For a counterclockwise rotation, braiding corresponds to the transformation 7 a — > +jb and 7b — > — 7 a while 
for a clockwise rotation we obtain the inverse transformation j a — > —7b and 76 — > +-f a 0- 




FIG. 3: Exchanging Majorana fermions in real space leads to a protected single topological qubit rotation termed braiding. In 
qubit space, this corresponds to a gate operation where the qubit states pick up a relative it/2 phase. 



For a system that supports a single MF per edge, the emergence of a unitary discrete symmetry O u with O 2 = I, 
can lead to an additional MF per edge. The two MFs per edge are labelled by the ±1 eigenstates of O u , leading to 
the following 4 edge MFs r y a ± and jb±. With the latter, we can define two zero-modes d± — (j a ± + ijb±)/V2 and 
two topological qubits with states {|1+) , |0+)} and , |0_)}. The accessible protected non-Abclian operations 

that we may perform within this four-fold degenerate Hilbert space are restricted by the simultaneous conservation 
of fermion parity and O u . Essentially, the possible operations are combinations of simultaneous or separate clockwise 
and counterclockwise braiding operations in each of the topological qubit spaces. Specifically, we have the following 
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four operations: 

A:{|1 + ),|0+),|1_),|0_)} H- {eW4| 1+))e -- /* | 0+ ) , |1_) , |0_» , (31) 

B:{|1 + ),|0 + ),|1_),|0_)} {|l + ),|0 + ),e-/ 4 |l_),e- iff / 4 |0_)}, (32) 

C : {|1+) , |0+) , |1_) , |0_)} -> {e"/ 4 |1 + ) , e""/ 4 |0 + ) , e i7r / 4 |1_) , e""/ 4 |0_)} , (33) 

D : {|1 + ) , |0 + ) , |1_) , |0_)} -+ {e"/ 4 |1 + ) , e^' A |0 + ) , e^ 4 |1_> , e+"/ 4 |0_)} , (34) 

presented in Fig. @] The operations A and £? correspond to braiding operations effected only on the + or the — 
topological qubits. These are single qubit operations. In contrast, if we effect braiding simultaneously in both qubits 
spaces we have two options. Either the direction of braiding is the same or opposite. These two possible topologically 
protected operations in the joint qubit space are described by the C and D configurations. The above set of protected 
operations do not suffice for performing universal TQC due to the Ising nature of the MFs [f| [5(|. Nevertheless, the 
presence of the additional topological qubit on the same wire can be useful for performing braiding operations. So 
far, several methods for performing braiding have been proposed, including networks of topological wires [H, [5l[ where 
neighbouring MFs can be controllably coupled in order to perform a MF exchange. In the present case, the additional 
O u protected MFs can constitute a reservoir of MFs that could reduce the number of complementary wires that one 
needs for performing adiabatic operations using these protocols. In addition, the presence of the extra pair of MFs 
can be also prominent for creating phase gate operations. A standard theoretical proposal for implementing a 
phase gate for two separated MFs, prescribes to bring the MFs to a finite distance in order to let them hybridize into 
a finite energy fermionic state. Due to the time evolution of the finite energy state, a phase gate operation will be 
implemented on the topological qubit when the MFs reseparate. In the presence of a hidden symmetry O u , one does 
not have to change the distance of the MFs any more. By controllably switching off the hidden symmetry O u , one 
hybridizes the two MFs of the same edge, for instance r y a ±, so to end up with a single MF. Depending on the details 
of the hidden symmetry breaking and restoration procedures, one may retrieve a phase gate operation. Of course a 
detailed investigation of these possibiltics is required. 

The alternative TQC routes described above depend delicately and crucially on the robustness of this hidden 
symmetry. As we have already mentioned, it is desirable to find a system that has a hidden symmetry related to a 
degree of freedom such as a band index. For example, in the case of a two-band topological superconductor where only 
intraband matrix elements appear in the Hamiltonian, the system splits into two irreducible subsystems similarly to 
the situation described above. As a matter of fact, multiband systems such as the Fe-bascd high-Tc superconductors 
[52| . offer a promising way out. The latter materials are supposed to exhibit intra-band superconductivity (usually 
a 4-band (53j | or a 5-orbital (lij picture is adequate) and consequently we may obtain a number of disconnected 
sub-systems. If we manage to render each of these superconducting sub-systems topological, we will be in a position 
to apply the topological quantum computing protocols discussed in the previous paragraph. Recently, a proposal 
concerning topological superconductivity based on iron-based superconductors has been put forward [55j . However, 
in that work an iron-based superconductor was used to induce superconductivity by proximity effects on a Rashba- 
semiconductor. Instead, the situation that I envisage involves intrinsic multiband topological superconductivity in 
the iron pnictide superconductor itself. 




FIG. 4: Allowed topologically protected operations for two topological qubits ± corresponding to the ±1 eigenvalues of a 
unitary hidden symmetry operator O u satisfying 0\ = I. The standard braiding operations A and B describe single qubit 
operations. C and D correspond to operations in the joint qubit space. In operation C, the ± pairs of Majorana fermions are 
braided in the same direction while in case D in the opposite. 
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VII. CONCLUSION 

I have performed a detailed analysis of the accessible topological superconducting phases that can occur from the 
combination of inhomogeneous Rashba spin-orbit coupling, magnetic order and superconductivity. By exploring the 
landscape of the possible topological phases I proposed new quasi-ld systems prominent for realizing MFs, based on 
Rashba spin-orbit coupling and T violating superconductivity, without the demand for any kind of magnetic order. 
Moreover, I pinpointed the significance of emergent unitary and anti-unitary hidden symmetries and revealed the 
topological implications that they lead to. Finally, I discussed alternative topological quantum computing pathways 
that open up in the presence of a unitary hidden symmetry and suggested that Fe-based multiband superconductors 
could be a potential candidate for these implementations. 

Note 

While the present manuscript was under completion, a preprint concer ning the classification of topological insulators 
and superconductors in the presence of a reflection symmetry appeared [561 ] . 
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